Transverse oscillations in solar spicules induced by 
propagating Alfvenic pulses 
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Abstract The excitation of Alfvenic waves in the solar 
spicules due to the localized Alfvenic pulse is investi- 
gated. A set of incompressible MHD equations in two 
dimensional x — z plane with steady flows and sheared 
magnetic fields is solved. Stratification due to gravity 
and transition region between chromosphere and corona 
are taken into account. An initially localized Alfvenic 
pulse launched below the transition region can pene- 
trate from transition region into the corona. We show 
that the period of transversal oscillations is in agree- 
ment with those observed in spicules. Moreover, it is 
found that the excited Alfvenic waves spread during 
propagation along the spicule length, and suffer efficient 
damping of the oscillations amplitude. The damping 
time of transverse oscillations elongated with decrease 
in kf, values. 
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1 Introduction 

Spicules are one of the most fundamental components 
of the solar chromosphere. They are seen in chro- 
mospheric spectral lines at the solar limb at speeds 
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of about 20 — 25 km s 1 prop agating from the chro- 
mosp here into the solar corona (jZaaarashvili fc Erdelvi 
20091) . Their diameter and length varies from spicule to 
spicule having the values from 400 km to 1500 km and 
from 5000 km to 9000 km, respectively. The typical life- 
time of them is 5 — 15 min. The typical electron density 
at heights where the spicules are observed is approxi- 
mately 3.5 x 10 16 — 2x 10 17 m~ 3 , and their temperatures 
are e stimated as 5000 
2000). 



8000 K (Beckers 



Kukhianidze et al. ( 20061) ; Zaqarashvili et al 



196 



Sterling 



(|2007j ) by analyzing the height series of Ha spectra 
in solar limb spicules observed their transverse oscilla- 
tions with the esti mated period of 20 — 55 and 75 — 110 
s. More recently, lEbadi et al.1 (|2012al ) based on Hin- 



ode /SOT observations estimated the oscillation period 
of spicule axis around 180 s. 

Despite the large body of theoretical and observational 
works devoted to the spicules, their ejection mecha- 
nism is not clear yet. In other words, observations 
with high spatial resolutions are needed to distinguish 
their origin. The Alfven waves are usual candidates 



corona ( 


De Pontieu et al. 


2007: 


Hollwesl 


( 1982|) showed that the 



Alfven waves may be nonlinearly coupled to fast magne- 
toacoustic sh o cks, w hich may lead to spicule formation. 



Cargill et al.l (| 1997T) performed the numerical simula- 
tions of the propagation of Alfvenic pulses in two di- 
mensional magnetic field geometries. They concluded 
that for an Alfvenic pulse the time at which differ- 
ent parts of the pulse emerge into the corona depends 
on the plasma density and magnetic field properties. 
Moreover, they discussed that this mechanism can in- 
terpret spicule ejection forced through the transition 
region. 



Kudoh fc Shibata ( 19991 ) used the randor 



linear Alfvenic pulses and concluded that the transi- 
tion region lifted u p to more than ^ 5000 km (i.e. the 
spicule produced). iDel Zanna et al.l (|2005l ) studied the 



2 



propagation and evolution of Alfvenic pulses through 
coronal arcades. Their results showed a strong spread- 
ing of the initially localized pulses due to variations 
of Alfven speed with height and efficient damping of 
oscillations. In addition, they concluded that the back- 
ground density and Alfven s peed should be considered 
as key ingredients in models. iMurawski fc Zaaarashvili 



:ey 

1 20101 ) studied the upward propagation of a velocity 



pulse launched initially below the transition region. 
The pulse quickly steepens to a shock which may form 
the spicules. They concluded that such a model may 
explain speed, width and heights of classical spicules. 
Rae fe Roberts! (1982) illustrated that impulsive mo- 
tions may propagate along an intense tube. Informa- 
tion (in the form of a wave front) propagates at the 
subsonic and subAlfvenic tube speed. As the wave 
front of an impulsively generated disturbance propa- 
gates along the tube, it trails behind it a wake, which 
oscillates at cut-off frequency. The significance of cut- 
off frequency depends upon what are the circumstances 
creating the motions. If disturbances are impulsively 
generated, then wave propagation occurs and the ex- 
istence of cut-off freque ncy is manifest in t h e crea tion 
of an oscillatory wake. lHasan fc Kalkofen dl999l) ex- 
amined the generation of waves in vertical flux tubes 
in the magnetic network through the impulsive excita- 
tion by granules. The buffering action of a granule on 
a flux tube at a certain level excites a pulse that travels 
away from the source region with tube speed. After the 
passage of the pulse, the atmosphere oscillates at the 
cut-off period o f the mode, with amplitude that s lowly 
decays in time. IZaqarashvili fc Skhirtladzd (|2008l ) con- 
cluded that photospheric granulation may excite trans- 
verse pulses in anchored vertical magnetic flux tubes. 
The pulses propagate upward along the tubes while os- 
cillating wake is formed behind the wave front in a strat- 
ified atmosphere. The pulses carry almost all the en- 
ergy of initial perturbations, while the energy in wake 
oscillations is much smaller. 

In the present work we perform two-dimensional simu- 
lations of MHD equations and show that the Alfvenic 
pulse which is launched below the transition region can 
excite Alfven waves and reach the height where spicules 
are observed. Section 2 gives the basic equations and 
theoretical model. In section 3 numerical results are 
presented and discussed, and a brief summary is fol- 
lowed in section 4. 



2 Theoretical modeling 

We keep effects of stratification due to gravity in 2D 
x-z plane in the presence of steady flow and shear field. 



Since the excitation of incompressible Alfven waves due 
to Alfvenic pulses is the main goal of this work, so 
the continuity and energy equations are not considered. 
The ideal MHD equations used are as follows: 



<9v 



-Vp + /)g + -(VxB)xB, (1) 



dB 

~dt 



V x (v x B), 



P 



pRT 
P 



(2) 



(3) 



Here, p is the plasma density; T is the plasma temper- 
ature in Kelvin; R is the universal gas constant, and 
p is the mean molecular weight. We assume that the 
spicules are highly dynamic with speeds that are signif- 
icant fractions of the Alfven speed. Perturbations are 
assumed independent of y, i.e.: 

v = v k + v y (x, z,t)j 

B = B 0x (x,z)i + B az (x,z)k + by(x,z,t)j, (4) 
and the equilibrium sheared magne tic field is two 



dimensional and divergence-free as (jDel Zanna et al 
20051 IMurawski fc Zaqarashvilill201oT ): 



B 0x (x,z) = B e k » z cos[k b (x - a)} 
B 0z (x,z) = -B e~ kbZ sm[k b (x - a)}, 



(5) 



with a as the spicule radius. Since the equilibrium mag- 
netic field is force-free, therefore the pressure gradient 
is balanced by the gravity force, which is assumed to 
be g=-g k via this equation: 



-Vpo(z) + p (z)g = 0, 
and the pressure in an equilibrium state is: 
dz' 

Po(z) =Po exp ' 



Zr A(z'). 

The density profile is in the form of: 



(6) 



(7) 



Po( z ) = ex P 



with 



A(z) 



To(z) 



RTpjz) 

pg 



dz' 



(8) 



(9) 
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where the temperature profile is taken here as a 
smoothed step function, i.e.: 



To(z) 



(1 - d t )tanh(- 



(10) 



here dt = T c h/T c with T c h denoting the chromosphere 
temperature at its lower part. The symbol T c corre- 
sponds to the temperature of the solar corona that is 
separated from the chromosphere by the transition re- 
gion, which has the width of z w — 200 km, and is lo- 
cated at the zt = 2000 km above the solar surface. We 
put T ch = 15 x 10 3 K and T c = 3 x 10 6 K. 
Figure Q] shows the equilibrium density and magnetic 
field with respect to x, and z (all quantities in these 
plots are non-dimensional). 

The linearized dimensionless MHD equations with 
these assumptions are: 



dvy 
dt 



dby 
dt 



v - 



8Vy 



Po(z) 



t-) / sdb y dby 
B 0x {x, z)-j£ + B 0z (x, z)-^ 



db 



y_ 

dz 



= B, 



O.r 



R 



{x,z)—± : 
dz 



(11) 



(12) 



where the densities, velocities, the magnetic field, 
time and space coordinates are normalized to po (the 
plasma density at dimensionless z = 6), V^o = 
Bo/ \J PoPOi Bq, t (the time scale of transit Alfven time 
defined as r = cl/Vao), a (spicule radios), respectively. 
Also the gravity acceleration is normalized to a/r 2 . 
Eqs.QTJ and [l2]should be solved under following initial 
conditions: 



v y (x, z, t = 0) = A v exp 



(a; - x a ) 2 + (z - z ) 2 



b y (x, z, t = 0) 



Ab sin(7ra;) sin(7rz) 



(13) 



where the initial magnetic perturbation amplitude is 
set to be Ab = 10~ 7 , and (xq, Zq) determine the center 
of an initial pulse where it has the maximum speed. 
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Fig. 1 (Color online) Dimensionless equilibrium density 
and total magnetic field plots with respect to dimensionless 
x, and z space coordinates (kb = 7r/3). 



3 Numerical results and discussion 

To solve the coupled Eqs. [TlJ and [12] numerically, the 
finite difference and the Fourth- Order Runge-Kutta 
methods are used to take the space and time deriva- 
tives, respectively. We set the number of mesh-grid 
points as 256 x 256. In addition, the time step is cho- 
sen as 0.001, and the system length in the x and z di- 
mensions (simulation box sizes) are set to be (0,2) and 



•1 



(0,8). The parameters in spicule environment are as fol- 
lows: a (spicule radios)=1000 km, uj = 0.3a = 300km 
(the width of Gaussian packet), L=8000 km (Spicule 
length), v = 25fcm/s, n e — 11.5 x 10 16 m~ 3 , B — 
1.2 x 10" 3 Tesla, T = 14 000 K, g = 272 ms~ 2 , 
R = 8300 m 2 s _1 A; _1 (universal gas constant), Vao = 
77.5 km/s, = 0.6, r = 13 s , p = 1.9 X 10~ 10 /.-.am" 3 , 
Po = 3.7 x 10~ 2 Nm~ 2 , /x = 4tt x 10~ 7 TmA' 1 , 
z r = 6000 fcm (reference height), z^, = 200 fcm, 



2000 km, A r , 



7.75 km/s, xn 



1. 



3C 



0.5 



(|Ebadi et al.ll2012bl libadi fc Hosseinpourll2013h . 

Figure [2] shows 3D plots of the perturbed velocity 
with respect to x, z at t = 13s, £ = 130s, and t = 185s 
for kb = 7r/3. The initial pulse amplitude is 7.75 km/s 
and is located in x = 1000 km, and z = 500 km. It 
shows that the initially localized Alfvenic pulse, which 
is launched below the transition region, propagates up- 
ward. We prefer to show the 3D plots of the perturbed 
velocity in the case of k^ = tt/3. Putting kb = 7r/3 in 
the equilibrium magnetic field leads to a pulse which 
propagates slowly regarding the total propagation dis- 
tance. So that, analyzing the temporal behavior of the 
pulse becomes easy. As seen, the initially localized 
pulse during propagation suffers significant spreading 
along the spicule length due to the variations of Alfven 
speed with height. This leads to an efficient damping 
of the oscillations amplitude (phase mixing). 

Similar to figure in figure |3] the same plots of the 
perturbed velocity with respect to x, z at t = 26s, 
t = 78s, and t = 130s for kb = ft/3 are presented, 
but the initial pulse is now located at x = 2500 km, 
and z = 500 km. Obviously, the similar pattern of 
upward propagation of the initially localized Alfvenic 
pulses repeated here. 

Figure 2] illustrates variations of transversal compo- 
nent of the perturbed velocity with time in x = 1900 
km, and z = 2500 km for different values of kb- The ini- 
tial pulse may excite the wake in stratified atmosphere 
which oscillates at cut-off frequency. In the top-left 
panel (kb = tt/3) of figure SI the oscillatory wake be- 
hind the initial pulse is presented. The period of the 
oscillatory wake corresponds to the period of transverse 



oscilla tions of spicules observed by IZaaarashvili et al 
(J2007J). The same pattern is appeared in other pan- 
els but with different periodicity values. An another 
interesting result, which can be seen in different pan- 
els of figure H] is that the damping time of oscillations 
is elongated towards the smaller kb values. The phase 
mixing and stratification due to gravity are the respon- 
sible mechanisms in damping of the oscil lations, and 
their efficiency reduce in small kb values (lEbadi et al 



2012bUEbadi k Hosseinpourll2013| ). We have not found 



any observational signature of this result in the litera- 
ture (thanks to De Pontieu, B. for his useful comment 
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Fig. 2 (Color online) The 3D plots of the transversal com- 
ponent of the perturbed velocity with respect to x, z in 
t = 13s, t = 130s, and t = 185s for kb = 7r/3. 
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Fig. 3 (Color online) The 3D plots of the transversal com- 
ponent of the perturbed velocity with respect to i, 2 in 
t — 26s, t = 78s, and t = 130s for kb = tt/3. 



about this point). It is, in fact, a very difficult mea- 
surement. 

4 Conclusion 

The initially localized Alfvenic pulse propagation in a 
medium with steady flows and sheared magnetic field 
is investigated. We take into account the stratification 
due to gravity and the transition region between chro- 
mosphere and corona. An initially localized pulse can 
excite the Alfvenic waves in the medium. It is shown 
that these waves can penetrate from the transition re- 
gion into the corona. The period of transverse oscil- 
lations (wake) that are induced in the medium due to 
the propagation of the Alfvenic waves are in agreement 
with those observed in spicules. Moreover, It is dis- 
cussed that the excited Alfvenic waves during propaga- 
tion along the spicule length suffer efficient damping of 
the oscillations amplitude. This can be related to the 
variations of Alfvcn speed with height (phase mixing) 
and gravitational stratification. It is interesting that 
the damping time of oscillations is elongated with the 
smaller kt, values. In other words, the phase mixing and 
stratification efficiency reduce in small fcj, values. 
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Fig. 4 (Color online) The time variations of transversal component of perturbed velocity in detection point (x = 1900 
km, and z = 2500 km) for different values of k^. 



